Entanglement in a molecular three-qubit system 
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We study the entanglement properties of a molecular three-qubit system described by the Heisen- 
berg spin Hamiltonian with anisotropic exchange interactions and including an external magnetic 
field. The system exhibits first order quantum phase transitions by tuning two parameters, x and y, 
of the Hamiltonian to specific values. The three-qubit chain is open ended so that there are two types 
of pairwise entanglement : nearest-neighbour (nn) and next-nearest-neighbour (nnn). We calculate 
the ground and thermal state concurrences, quantifying pairwise entanglement, as a function of the 
parameters x, y and the temperature T. The entanglement threshold and gap temperatures are also 
determined as a function of the anisotropy parameter x. The results obtained are of relevance in 
understanding the entanglement features of the recently engineered molecular Cr-jNi-Cv?'^ -Cr^Ni 
complex which serves as a three-qubit system at sufficiently low temperatures. 

PACS numbers: 03.67.Mn, 03.67.Bg, 03.65.Ud, 64.70.Tg, 75.10.Dg 



I. INTRODUCTION 

Entanglement is a unique feature of quantum mechan- 
ical systems with no classical analogue. In an entangled 
state, two or more quantum particles have joint prop- 
erties in the form of non-local correlations rather than 
individual identities. Entanglement is known to be a key 
resource in quantum information processing (QIP) tasks 
such as quantum computation, teleportation and cryp- 
tography Implementation of QIP protocols requires 
the assembly of multi-qubit systems with the potential 
for generating controlled entanglement. Natural exam- 
ples of qubits, which are two-level systems, include spin- 
i particles, photons with two states of polarization and 
trapped ions with two atomic states. In recent years, 
molecular nanomagnets have been proposed as appro- 
priate candidates for qubit encoding and manipulation 
[3, S|- A specific example is provided by antiferromag- 
netic (AFM) Cr-jNi rings which reduce to effective spin- 
\ systems at low temperatures. Each octagonal ring 
consists of one Ni"^^ and seven Cr^+ ions with AFM 
coupHng between neighbouring ions. A variety of exper- 
imental techniques have been used to characterize the 
rings. The rings have spin-i ground states and behave 
as qubits at sufficiently low temperatures as the excited- 
state multiplets remain unoccupied. Also the rings have 
been demonstrated to possess long decoherence times, an 
ideal requirement for several QIP tasks. 

Recently, Timco et al j3| have engineered a coherent 
coupling between two Cr^Ni rings, serving as molecular 
spin qubits, via a central Cv?^ ion which acts as a third 
qubit. The CrjNi-Cu'^^-CrTNi complex is equivalent 
to a three qubit system with the Cu^"*" ion serving as 
a "linker". The coupling between the spins of the rings 




Figure 1: A molecular three-qubit system in which the 
qubits A and B represent two Cr-jNi rings (see figure 

1(a) of reference Q) and qubit C represents the 
bridging ion Cv?'^ . The rings and the ion are effective 
spin-i systems represented by solid arrows. 



is tunable by a proper choice of the linker. In a mi- 
croscopic approach, the spin Hamiltonian describing the 
three-qubit system can be written as 01 



H = 



H 



AC 



H 



BC 



(1) 



where the labels A, B and C correspond to the two rings 
and the magnetic Hnker respectively. The terms H"^ and 
individually describe the CrjNi rings: 
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with z along the ring axis. The successive terms in the 
Hamiltonian correspond to isotropic exchange (Ji), axial 
crystal field (di), dipole-dipole coupHngs {Hdip) between 

eight individual spins S i and the Zeeman coupling to the 
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magnetic field B with gi being the gyromagnetic tensor. 
The term H'-' in equation (1) is 



H'" — B^cu - S c 



whereas the terms H^^ /H^^ are: 



H 



AC _ ttBC _ Tf 



' Cr + S Ni 



(3) 



(4) 



where the spins S cr and S Ni correspond to 5 1 and 
S's in their respective rings as these spins are located 
on the edge of the octagon bound to the Cu link. Since 
J' « Ji^s, the intra-ring exchange constants, the low- 
temperature behaviour of the CrT N i-Cv?^ -Cr-j N i com- 
plex is determined by the splitting of the lowest eight 
energy levels. The behaviour can be reproduced in terms 
of an effective three-spin Hamiltonian [Sj: 

= J ^ Si.Sc + ^j.bB . ^ ^-Si 

i=A,B i=A,B,C 

+Dex ^ C^Si^zScz — Si^xSc,x — Si^yScy) (5) 
i=A,B 



where S a,b,c represent spin-i operators, J is the 
strength of the effective Cw-ring isotropic exchange, 

~g A,B are the 5-tensors of the ring ground doublet, 

9c— 9 Cu, and Dex is an effective Cu-ring axial ex- 
change originating from the projection of the rings' dipo- 
lar and crystal-field anisotropies. Equation (5) represents 
the CrrNi-Cu'^^-CrTNi system as a linear chain of three 
coupled qubits with open boundaries. The three-qubit 
system has ground and thermal states which are entan- 
gled. One can focus on two types of entanglement : pair- 
wise, i.e. between two qubits and three-party entangle- 
ment involving all the three qubits. The Greenberger- 
Horne-Zeilinger (GHZ) and Werner (W) states 0, S 01 
defined as 



|Gi/z> - ^ (I TTT) + Uii)) 
TTi> + l TiT> + UTT» 



(6) 



represent two fundamentally non-equivalent entangled 
states of three qubits. In the first case, the pairwise en- 
tanglement for all the qubit pairs is zero and one has gen- 
uine three party entanglement known as the residual en- 
tanglement. The nomenclature arises from the Coffman- 
Kundu-Wootters (CKW) inequality Q for a three qubit 
system given by. 



Tl > T2 



(7) 



where ti represents the one-tangle corresponding to the 
entanglement between the ith qubit and the rest of the 
system and Cfj is the square of concurrence, a measure 
of the entanglement between the ith and jth qubits. The 
one-tangle ri is determined as ti = Adetp'^^^ where p'^'' 



is the single-site reduced density matrix. The residual 
entanglement is given by the difference between ti and 
T2 and hence provides a measure of quantum correlations 
which cannot be expressed in terms of pairwise correla- 
tions. The GHZ state has the maximum possible value 
of 1 for the three-party (residual) entanglement. The W 
state, on the other hand, possesses only pairwise entan- 
glement between all qubit pairs and the magnitude of the 
residual entanglement is zero. Timco et al Q have pro- 
vided a prescription for the generation of GHZ and W 
states using a sequence of microwave pulses applied to 
the molecular three-qubit system. 

In this paper, we study the entanglement properties 
of the ground and thermal states of the molecular three- 
qubit system described by the reduced Hamiltonian in 
equation (5) . We specially focus on the variation of en- 
tanglement measures as a function of the parameters of 
the Hamiltonian. Wang et al (1] have earlier studied 
thermal entanglement in the three-qubit Heisenberg XXZ 
model. The Hamiltonian considered by them satisfies 
periodic boundary condition and includes anisotropic ex- 
change interaction and magnetic field terms. The molec- 
ular three qubit system considered in this paper has the 
structure of an open chain and the entanglement features 
turn out to be different from those of the three-qubit 
Heisenberg ring. The experimental demonstration that 
the coupling between the molecular spin clusters can be 
controlled without disturbing the intra-cluster interac- 
tions provides the impetus for characterizing the entan- 
glement properties of the molecular three-qubit system. 



II. GROUND STATE ENTANGLEMENT 

We consider the molecular three-qubit system to be in 
an external magnetic field pointed in the z direction. The 
Hamiltonian (equation (5)) then reduces to 



H = J ^ S,.Sc+911bB ^ St 

i=A,B i=A,B,C 



+D. 



ex 

i=A,B 



{2Si^zSc,z — Si^xSc,x — Si^yScy) (8) 



This can be rewritten in the form. 



H^H/J={1 + 2x) {Sa^.Sc^z + Sb^.Sc,.) 
1 

+ 2 

+y{SA^z + SB^z + Sc,z) (9) 



(1 x) (yS\S(-, + Sj^S'^ + SgSfj + SgS'^) 



where x = Dex/ J,y = g^^BB/J and S'+, S~ are the rais- 
ing and lowering operators. Since the ^-component of the 
total spin, 5"*°*, is a conserved quantity, the eigenvalue 
problem can be solved in the separate subspaces corre- 
sponding to the different values of 5**°* . The eigenvalues 
and the eigenstates are given by. 



S\ 



tot _ _|_3 
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IV'4) = ni)-s{x)\ nt) 

Ei^\{2y~U^{x)} 



iTT» 



(13) 



5*' 



m = Au^i\m)-BMin) + \ni)) 

E, = ^{~2y-U+{x)} 



1^6) = (I uT)-5(a:)im) + ini)) 



(14) 



(15) 



|V'7) = ;^(-|UT> 



Tii» 



(16) 



5"*°* = - 



1^8) = I III) 

Es^^{l + 2x- 3y) 
In the above equations, 

U±ix) = I 1 + 2a; ± v^S {Ax^ - 4a: + 3)| 



Rix) 



S{x) 



2(-l + a;) 
2(-l + a;) 



(17) 

(18) 
(19) 

(20) 



A{x) 



B{x) 



2(-l + a;) 



-U-{^) 
2(-1 + .t) 



(21) 



(22) 



We first consider the case of zero magnetic field (y=0). 
The eigenvalues then become 



El — Es 



l+2x 



^8 — 2 

E2=Et = Q 



E,=E, = ^^ 
E, = E, = 



(23) 



We assume x to range over both positive and negative val- 
ues. Figure 2 shows a plot of the two lowest energy levels 
of the Hamiltonian (equation (9)) versus the parameter 
X. Each energy level is doubly-degenerate. The nature 
of the ground states change at x = — 2, bringing about a 
first-order quantum phase transition (QPT). When x is 
< — 2, the ground states are the separable states jV'i) and 
IV's)- When a; is > —2, the doubly-degenerate ground 
state is described by the wave functions |?/'3) and IV's)- 
At a: = 1, however, the Hamiltonian (9) becomes Ising- 
like, i.e., loses its quantum character and the degenerate 
ground states, | tit) and | iti), are separable. We now 
discuss the entanglement properties of the ground states. 
Because of the degeneracy, the ground state density ma- 
trix describes a mixed state with 



IV'5)(V'5|) 



(24) 



The reduced density matrix p^, — A,B,C) is ob- 
tained from p by tracing out the spin degrees of freedom 
associated with the spins which are not located at the 
sites i and j. The reduced density matrix in the stan- 
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Figure 3: Variation of concurrences (a) Cac and (b) 
Cab versus x for y — 0. 



Figure 4: Plot of the two lowest energies (i)i?8 and 
{ii)E5 versus x for y — 0.5. 
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Figure 5: Plot of the lowest energy levels {i)Es, (ii)i?5, 
(iii)i?3 and (iv)i?i as a function of y for a: = 0.5 



dard basis, {| tt), I Ti), I it), I ii)}, has the structure 



/ u 

wi y* 

y W2 

VO V, 



(25) 



The concurrence Cy , a measure of the entanglement be- 
tween a pair of spins at sites i and j, is given by [gl. flo|. 



Cij — 2 max 



(o,|y| 



uv) 



(26) 



Figures 3(a) and 3(b) show the variation of Cac and Cab 
versus x. The analytical expressions for the concurrences 
are: 



Cac = Cbc = 2 max (O, | ^ | - 
Cab — 2 max (^0, 



1 

I A' 



IT- 



(27) 



The variation of Cbc as a function of x is identical with 
that of Cac- We remind ourselves that A and B are 
the boundary spins and C the central spin. A jump in 
the magnitude of the concurrence indicates a first order 
QPT [ni, [3, which, as already mentioned, occurs 
at a; = —2. Cac and Cbc both become zero at a; = 1 due 
to the separability of the ground state density matrix and 



then rises as x is increased to attain a saturation value 
Cac = g^^yl for large x. The entanglement between 
the boundary spins, however, has a non-zero value only 
for negative values of x and that too in a restricted range 
of X values. 

We next consider the case of non-zero magnetic field 
{y ^ 0). Figure 4 shows the plots of the two low- 
est energies, E5 and Es, versus a; for y = 0.5. One 
finds that a first-order QPT occurs at a specific value 

of a; = aic ^= — — indicating a change in the 

nature of the ground state. Figure 5 shows the variation 
of Ec, and Eg as a function of y for x = 0.5. Again, one 
notes the occurrence of a first-order QPT at a specific 
value of y = yc- The external magnetic field removes 
the ground state degeneracy of the zero-field case and 
the three-qubit system has a unique ground state. Fig- 
ures 6(a) and 6(b) show the variation of the concurrences 
Cac and Cab versus x for y — 0.5. In this case, the an- 
alytical expressions for the concurrences are; 



Cac ^ Cbc = 2 max {0,\^ 
Cab = 2 max (O, |^|) 



(28) 



In zero magnetic field (y = 0), the next-nearest- 
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Figure 6: Variation of concurrences (a) Cac and (b) ^ig^'^e 7: Variation of concurrences (a) Cac and (b) 

Cab versus x for y = 0.5. versus y for a: = 0.5. 



neighbour (nnn) concurrence Gab has non-zero values 
only in a restricted range of negative x values (figure 
3(b)) whereas in the presence of a magnetic field (y 7^ 0), 
Cab is non-zero in a range of both negative and positive 
X values. The magnitude of the nnn entanglement is less 
than that of the nn entanglement for both y — Q and 
y ^ 0. As y increases, one finds that Xc, the first or- 
der QPT point shifts towards more positive values. For 
sufficiently high values of y, entanglement exists only for 
positive values of x. This is so provided y is less than 
the critical value yc (which depends upon x) at which a 
first-order QPT takes place to a separable ground state. 
Figures 7(a) and 7(b) show the plots of the nn and nnn 
concurrences, Cac{— Cbc) and Cab respectively versus 
y ior X — 0.5. The concurrences have constant values 
for y < yc, the QPT point, and jump discontinuously 
to zero values at y = yc- As x increases, the value of 
yc also increases. In the case of non-zero magnetic field, 
y ^ 0, the ground state is non-degenerate and the den- 
sity matrix represents a pure state. In this case, the 
one-tangle ti (defined in equation (7)) can be calculated. 
For any choice of the central spin, the residual entangle- 
ment involving three spins is found to be zero so that 
only pairwise entanglement exists in the ground state. 
The ground states thus belong to the class of W rather 



than GHZ states. 



III. THERMAL STATE ENTANGLEMENT 

We next discuss the finite-temperature entanglement 
properties of the molecular three-qubit system. The ther- 
mal density matrix p{T) = exp {—(]H){(3 = j^,kB = 1) 
now replaces the ground state density-matrix with Z de- 
noting the partition function of the system. The reduced 
density matrix Pij(T) has the same form as in equation 
(26) with Cij (T) given by [Hi 

a, (T) - I max (0, \y{T)\- ^u{T)v{T)) (29) 

For the three-qubit system, the thermal density matrix 
is 

1 ^ 

/'(^)^^E^'^P("'^^fc)l^fc)<^'^l (30) 
fe=i 

where the |i/'fc)'s and iJ^'s are given in Eqs. (10)-(17) for 
the general case y 0. When y = 0, the energy eigenval- 
ues are as shown in equation (23). We first consider the 



6 



case of zero magnetic field (y = 0) . The matrix elements 
w, V and y of the reduced density matrix pAc{T) are 



1 -E3 
e ^ 



2R 



1 

52 



-E4 



25 



(31) 



(32) 



The reduced density matrix pbc{T) has the same ma- 
trix elements as in equations (31) and (32). For the nnn 
concurrence, Cab, the matrix elements of the reduced 
density matrix are : 



i?2 



y 



52 ^ 

p 1" 



— ?re ^ — 1 



(33) 



(34) 



Figures 8(a) and 8(b) show the plots of Cac and Cab 
respectively as a function of x for different values of the 
temperature T. As T increases, the range of x values 
for which Cac 7^ shifts towards more positive values. 
Figures 9(a) and 9(b) show Cac and Cab versus T for 
negative values of x. One can obtain similar plots for 
Cac when x is > 0. For both the nn and nnn entangle- 
ments, one can define threshold temperatures T^'' and 
T^^ respectively, beyond which the concurrences [1, 

have zero values. Figures 10(a) and 10(b) show how T^'' 
(2) 

and vary with x for different values of y. 

For non-zero magnetic field, y ^ 0, the matrix elements 
of the reduced density matrix pAc{— Pbc) are given by. 



El 1 £2 1 _iE3 1 
T 4- — p T -4- p T 4- f 

^2 ^ ^ B^ 



(35) 




Figure 8: Variation of concurrences (a) Cac and (b) 
Cab versus x at different temperatures T with y — 0. 
The different temperature values are (a) (i) T = 0, (ii) 
T = 0.5, (iii) T = 1.5 and (b) (i) T = 0, (ii) T = 0.1, 
(iii) T = 0.3 



Figures 11(a) and 11(b) show the plots of the nn and nnn 



/ Eg 1 _ Et 1 _ E5 1 Eg \ 

g ~ -) e ~ -] e ~ H e ~ (36) 

\ 2 A"^ B"^ J entanglements, Cac and Cab respectively, versus a: for 

y — 0.5 and at different values of T. Figure 12(a) shows 

_ ^ the plot of Cac versus T for different negative values of x 

^ s ^ '^^^ R S ^6 A (37) ■^ith y = 0.5. Similar plots are obtained in other ranges 

A^^ ^ ~ ^ ~ 1^^ ^ ~ ^ ) of X values. Figure 12(b) shows how Cab varies as a 

function of T. 

The corresponding matrix elements for the nnn reduced 

density matrix are: 1 1 1 ^ ^.i ^ 1 

■' We lastly calculate the entanglement gap temperature 

^2 \ [i^^ [1^^ Te, as a function of x for both zero and non- 

^ ^ ^ ^ j zero y. Te is determined from the relation ?7(Tb) = E'sep, 

where U{T) ^= —-^^^ is the thermal energy at temper- 
^2 ^2 \ ature T and Esep is the ground state energy of the clas- 

^ 52 ^ ^ j ^^^^ sical spin model corresponding to the three-qubit Hamil- 

tonian in equation (5). Esep can be easily calculated, 
e.g., Esep — \ (—1 — 2x — y) for x > 0. For temperature 
T < Te, the thermal state is entangled. Figure 13 ex- 



y = y =-2\^ T+e 



1 



^ Eg \ (40) hibits the variations of Te versus x for different values of 
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(b) 



Figure 9: Variation of concurrences (a) Cac and (b) 
Cab versus T for different negative values of x with 
y = 0. The different x values are (i) x — —1.5, (ii) 
X = —1, and (iii) x — —0.55 for both (a) and (b). 



1.5 - 
1.0 - 
0.5 

0.0: 

0.4 
0.3 
0.2 
0.1 
0.0 



(a) 













\ (iii) 
Vii) \ 















(b) 

Figure 10: The threshold entanglement temperature (a) 



T^^-* versus x for nn entanglement and (b) T^' versus x 
for nnn entanglement with (i) y = 0.1 (ii) y = 0.5 and 
(iii) y — I- Tj^"* has a very weak dependence on the 
values of y. 



(2) 



IV. SUMMARY AND DISCUSSION 

In this paper, we have obtained quantitative measures 
of pairwise entanglement in a molecular three-qubit sys- 
tem as a function of two parameters x and y. The system 
represents the recently engineered Cr'jNi-Cu'^^-Cr'jNi 
complex consisting of two Cr^Ni rings coupled via a cen- 
tral Cv?'^ ion. The parameters x and y appearing in the 
qubit Hamiltonian (equation (9)) have their origins in 
an effective Cw-ring axial exchange due to the projection 
of the rings' dipolar and crystal-field anisotropics and 
an external magnetic field respectively. Timco et al Q 
have provided an experimental demonstration that the 
coupling between molecular spin clusters can be manipu- 
lated by altering the nature of the Hnker ions. This opens 
up the possibility of chemically controlling the generation 
of entanglement in spin systems. The molecular three- 
qubit system studied in this paper belongs to a family of 
clusters with AFM exchange interactions between the nn 
ions and a spin-i ground state. The simplest case is that 
of a finite chain with an odd number of S" = ^ spins and 
dominant AFM interactions between the nn spins. An 
alternative way of obtaining an S = ^ ground state is 



to replace a single spin in an AFM chain, containing an 
even number of spins, by a spin of different magnitude 
such that the ground state spin is of magnitude ^ . The 
Cr^Ni ring provides an example of the latter possibility. 

In general, the arrangement of spins in a chain can be 
either Hnear, or cyclic. The molecular three-qubit system 
studied in this paper is a linear-chain complex, whereas 
the three qubit chain studied in Q is cyclic in nature. 
Cyclic spin chains with an odd number of antiferromag- 
netically coupled spins have degenerate ground states due 
to magnetic frustration. For a three-spin cyclic chain, 
the anisotropic Heisenberg XXZ model has a four-fold 
degenerate ground state in zero magnetic field fS] . In the 
presence of an external magnetic field, the ground state 
is doubly degenerate. The effective three-qubit Hamil- 
tonian (equation (9)) has the form of the anisotropic 
Heisenberg XXZ Hamiltonian but with a linear i.e. an 
open-ended structure. In this case, the ground state is 
doubly degenerate in zero magnetic field {y = 0) and non- 
degenerate when y ^ Q. The cyclic chain has a greater 
ground state degeneracy because of frustration. 

There are prominent differences in the entanglement 
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Figure 11: Variation of concurrences (a) Cac and (b) 
Cab versus x at different temperatures T with y = 0.5. 
The different temperature values are (a) (i) T = 0, (ii) 
T = 0.5, (iii) T = 2 and (b) (i) T = 0, (ii) T = 0.5 
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(b) 

Figure 12: Variation of concurrences (a) Cac and (b) 
Cab versus T for different negative values of x with 
y = 0.5. The different values of x are (i) x = —1.5, (ii) 
X = —1, and (iii) x — —0.5 for both (a) and (b). 



features of cyclic and linear spin chains. As shown in 
in the case of the AFM cyclic XXZ model, there is 
no pairwise entanglement, as measured by concurrence, 
for all values of the anisotropy constant. At T = also, 
the concurrence is zero for the AFM case. In contrast, 
the AFM linear chain has pairwise entanglement at both 
T = and T (figures 3 and 8). One now distinguishes 
between nn and nnn entanglements. In the AFM case, 
the nnn concurrence Cab is zero at zero and finite tem- 
peratures whereas the nn concurrence Cac{ — Cbc^ is 
non-zero at both T — and T 7^ 0. Comparing figures 3 
and 6, one finds that on inclusion of the magnetic field the 
range of x values for which Cab (the nnn concurrence) 
is ^ is considerably extended. For a specific value of 
X, there is, however, a critical value of yc of y such that 
the pairwise entanglement vanishes when y > yc- Earlier 
studies [1, fioj have shown that the entanglement between 
two spins in an AFM chain can be increased by raising 
the temperature or the external magnetic field in specific 
ranges. This is true for our three-qubit system also. In 
figures 12(a) and 12(b), the curve (i) shows the increase 
of both Cac and Cab with temperature T. We have 
further shown that only pairwise entanglement exists in 



the ground state with y 7^ 0, i.e., there is no three-qubit 
entanglement as exists in the GHZ state (equation (6)). 
One interesting feature of the linear three-spin chain re- 
lates to the variation of the threshold entanglement tem- 
peratures T^^ and T^^ versus x for different values of 

y. As shown in figure 10(a), the T^^ versus x plot de- 
pends weakly on the values of y. The threshold temper- 
ature Tq \ for nnn entanglement, however, varies more 
prominently with y. In the case of the cycHc chain, the 
single threshold temperature depends on both x and y. 
As shown in figure 13, the plots of the entanglement gap 
temperature, Te, versus x are different for different val- 
ues of y. In fact, Te has a non-monotonic dependence on 
the values oiy (the y = I curve lies in between the y = 0.1 
and y = 0.5 curves). One further notes, from figures 10 
and 13, that the entanglement gap and threshold temper- 
atures are different for the same values of the parameters 
X and y. In fact, one finds that T^^'' < Te < T^\ Fig- 
ures 4-7 and Figures 11-12 have been obtained by fixing 
either a; or y at a specific value. The observations are, 
however, general in nature and hold true in extended 
ranges of x and y values. In the model studied by us, we 
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0.0 0.5 1.0 1.5 2.0 



X 

Figure 13: Variation of the entanglement gap 
temperature, Te , versus x for different magnetic fields. 
The different values of y are (i) y = 0, (ii) y = 0.5, and 
(iii) y = 1 



have assumed that the gyromagnetic factors gA, 9b, 9c 
are of equal amplitude 9 (equation (8)). In the case of 
the engineered three-qubit system, the diagonal tensors 
9A,B and 9c are different. Assuming 9a =_9b 7^ 9c, 
{9c{zz) = 2.07, 9a,b{zz) = 1-79, as quoted in 01), we find 
no qualitative changes in the results reported in sections 
2 and 3. It will be of interest to study the general case 
of the magnetic field pointing in an arbitrary direction. 
The three-qubit molecular cluster exhibits first-order 



QPTs at specific values of x and y. In figure 3(a), the 
QPT at a; = — 2 separates two phases, for x < —2 the 
ground state has no entanglement whereas for — 2 < a; < 
1, the ground state, described by the mixed state in equa- 
tion (24), has pairwise entanglement. Similarly, as shown 
in figure 7, a first order QPT occurs at y = j/c- The 
point a; = 1 is of special interest as the ground and ther- 
mal states become separable at this point. The threshold 
entanglement temperatures, T^^ and drop sharply 
to zero at a; = 1. The first-order transition points can 
be shifted by changing the parameters x and y. For ex- 
ample, the transition point Xc can be shifted towards 
higher values by increasing y. The first order QPTs 
are marked by discontinuities in the magnitude of both 
the nn and nnn concurrences associated with the ground 
states. The molecular three-qubit system, CttNi-Cv?^- 
CrrNi, has been specifically engineered with QIP appli- 
cations in mind. Since entanglement is a fundamental 
resource in such applications, a knowledge of its depen- 
dence on the relevant parameters of the system will be 
of use in the designing and implementation of QIP pro- 
tocols. With possibilities for controUing the couplings 
in molecular qubit systems 0| and realizations of spin 
Hamiltonians in optical lattices jl^ , some of the theoret- 
ical results could be observed in actual experiments. 
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